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ABSTRACT 

The Submillimetre Common User Bolometer Array (SCUBA) instrument has been op- 
erating on the James Clerk Maxwell Telescope (JCMT) since 1997. The data archive 
is now sufficiently large that it can be used to investigate instrumental properties and 
the variability of astronomical sources. This paper describes the automated calibra- 
tion and reduction scheme used to process the archive data with particular emphasis 
on 'jiggle-map' observations of compact sources. We demonstrate the validity of our 
automated approach at both 850- and 450- urn and apply it to several of the JCMT sec- 
ondary flux calibrators. We determine light curves for the variable sources IRC+10216 
and OH 231.8. This automation is made possible by using the ORAC-DR data reduc- 
tion pipeline, a flexible and extensible data reduction pipeline that is used on UKIRT 
and the JCMT. 

Key words: methods: data analysis - astronomical data bases: miscellaneous - in- 
frared: stars - telescopes - ISM: individual(IRC+10216) 



1 INTRODUCTION 

The Sub millimetre Common User Bolometer Array 



(SCUBA) (Holland et al. 1999) consists of two arrays of 
bolometers (or pixels); the Long Wave (LW) array has 37 
pixels operating in the 750- and 850- u.m atmospheric trans- 
mission windows, while the Short Wave (SW) array has 91 
pixels for observations at 350- and 450- u.m. Each of the pix- 
els has diffraction-limited resolution on the telescope (ap- 
proximately 14.5 and 8 arcsec FWHM respectively), and 
are arranged in a close-packed hexagon. Both arrays have 
approximately the same field-of-view on the sky (diameter 



not feasible to reduce the data manually and an automated 
approach must be used. 

There are three reasons why automated use of the 
archive is required. Firstly, it provides a unique resource for 
the investigation of instrument and telescope parameters, 
leading ultimately to an optimisation of observing practices, 
secondly the proposed Virtual O bservatories will requ ire re- 
duced and calibrated data (e.g. Jenness et al. 2002), and 
thirdly it allows monitoring of variable sources, especially 
extragalactic radio sources that are observed frequently for 
pointing purposes. In this latter respect the prese nt paper 



of 2.3 aieiiiiu), and can be used simulLaiieuusly by means uf 



a dichroie 



forms a companion to ou r recent work on this topic ( Robson 



Stevens & Jenness 2001) 



SCUBA was commissioned in late 1996 and there now 
exists a large searchable archive of observations^]. As of 
February 2001 the archive contained approximately 80,000 
separate observations. These data provide an excellent his- 
torical record of the telescope and instrument behaviour as 
well as being an archive of astronomical data. The major 
challenge of an archive of this size is the data processing. In 
order to extract useable historical data from the archive it is 



In order to determine its accuracy, any automated tech- 
nique under consideration must be used to process a stan- 
dard calibration source of known properties. We present our 
choice of such a reduction method below, using Uranus as 
the primary calibrator and then apply the technique to Mars 
and several compact secondary calibration sources as a con- 
sistency check. 



The public archive is hosted by the Canadian Astronomy Data 
Center 
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2 OBSERVATIONS 

In order to test any automated reduction system it is 
important that a suitable dataset is extracted from the 
archive. This dataset should contain frequent observations 
of a bright, compact source of known flux taken with a 
standard observing mode. Although Mars is the primary 
calibrator at submillimetre wavelengths, Uranus is a bet- 
ter choice for our approach because of its compactness rel- 
ative to the JCMT beam (semi-diameter approximately 1.7 
arcsec over the period in question). The best sampled , ho- 
mogen eous data set comes from jiggle-map mode (e.g., Hol- 
land et al. 1999) observations of Uranus that are performed 



at 850/450-u.m for calibration purposes, and at 850-u.m for 
pointing checks. We have searched the archive for all jiggle- 
map observations of Uranus; between April 1997 (the date 
when the SCUBA data headers were finalised) and May 2001 
there were approximately 1500 850- \im observations that 
were performed in conditions suitable for accurate extinc- 
tion correction, with 300 of these including data at 450-u.m. 
These data use chop throws of between 60 (all the pointing 
observations) and 120 arcsec (most of the maps). In Octo- 
ber 1999 both the 450- and 850- \im filters were updated to 
new wideband versions (denoted with a 'W suffix cf. 'N' 
for the previous filters). We include analysis for both filter 
combinations. 



3 DATA PROCESSING TECHNIQUES 

All SCUBA data are reduced using the SURF (SCUBA User 



Reduction Facility) pack age (Jenness & Lightfoot 1998; Jen 



ss & Lightfoot 2000 ) in conju nction with the Starlink 
software environment (Bly 2000). To allow us to process 



the many thousands of observations described in this paper 
we automated the reducti on using the ORAC data reduc- 
tion pipeline (OR AC-DR) (Economou et al. 1999; Jenness & 
Economou 1999| ). The general reduction method for SCUBA 
data is documented elsewhere (e.g., |5andell et al. 200l| ). In 
this section we comment briefly on the steps that we have 
tailored to our automated procedure, or which require some 
discussion in this context. 



3.1 Extinction correction 

Once the data have been demodulated, corrected for beam 
switching and flatfielded a critical calibration step is to cor- 
rect for the atmostpheric attenuation. The difficulties of 
calculating the atmospheric opacit y at submillimetre wave - 
lengths are well documented (e.g., |5tevens & L jtoteon_ji)94|} 



Archibald et al 



but re cent progress has been presented by 
(2002) who use skydips taken at both the JCMT and the 
Caltech Submillimeter Observatory (CSO) to provide a reli- 
able record of the opac ity spanning a period o f several years. 
We refer the reader to Archibald et al. (2002) for a detailed 



discussion but give a brief outline below. 

Automated reduction relies critically on the quality of 
the opacity measurement - simply applying the nearest sky- 
dip or CSO tau measurement to the data is dangerous for a 
number of reasons. Specifically, SCUBA skydips are usually 
not taken frequently enough to be useful in variable condi- 
tions. In addition, although the CSO tau measurements are 



taken much more frequently (every 10 mins), the instrumen- 
tal noise on the CSO tau meter is quite high, it produces er- 
roneously large readings occasionally and operates at a fixed 
azimuth. 

With these points in mind we p roceeded as follows: a 
polynomial was fitted to the re-scaled ( Archibald et al. 2002 ) 
CSO tau data in regions where the SCUBA skydips were 
in reasonable agreement; times when the CSO tau meter 
either did not agree with the SCUBA value or when the 
signal was varying rapidly were not included in these fits. 
This method ensured that the atmosphere seen by the tau 
meter was the same as that seen by SCUBA. In practice, 
the pipeline ignores data taken in periods when no fitted 
extinction value is available. Whilst this means that no data 
are reduced when, for example, the CSO tau meter is broken 
it is a necessary requirement of the automation process. The 
error of the 225 GHz opacity from the fit is usually of order 
0.005; if this error is assumed to be the dominant error in the 
extinction then for a source at airmass 1.5 this contributes 
to the final result an error of approximately 3 per cent at 
850- \im and 20 per cent at 450- vim. 



3.2 Despiking 

Occasionally, the time-series data contains spikes (for ex- 
ample from cosmic rays) and these must be removed. In the 
simplest case this involves analysing each bolometer in turn, 
removing data that are significantly different to the mean. 
The complication for SCUBA mapping is that each bolome- 
ter sees different parts of the target during the observation 
( "jiggling" on and off a source) such that a simple mean no 
longer provides valid statistics. For long observations it is 
possible to compare data points that are from similar points 
on the sky (effectively equivalent to subtracting source sig- 
nal from the data prior to despiking) but for short observa- 
tions, such as those from bright calibrators described below, 
there are simply too few data points for a reliable detection 
of spikes when near the source in question. For this work, 
we have compromised by performing the despiking for each 
bolometer independently but with a large enough clipping 
level (5 sigma) that we can detect the largest spikes with- 
out fear of clipping real source structure. This approach does 
mean that we are still susceptible to small spikes in the data 
but the overall effect on the integrated sum will be minimal 
when averaged over all the observations in this sample. 



3.3 Sky-noise removal 

Nodding and chopping do not remove all the atmospheric 
emission noise from the data so an additional step is re- 
quired that takes into account the fact that the source sig- 
nal is constant in time and the sky-noise is varying in time 
but coherent across the array. Removal of sky-noise from 
jiggle maps requires knowledge of regions of the map that 
are uncontaminated by source flux. Our criterion of selecting 
only compact sources for analysis simplifies this procedure 
enormously. For pointing observations it is safe to use the 
median of all the bolometers to determine the sky-noise sig- 
nal; the pointings use an azimuthal 60 arcsec chop-throw 
so that each map contains the source at the centre of the 
map and two negative images with half the flux. The latter 
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images are present because the chop throw is smaller than 
the field-of-view with the result that the negative beams of 
some of the 'sky' bolometers fall onto the central source. 
For standard mapping observations the chop throw is gen- 
erally chosen to be larger than the field-of-view in order to 
avoid this effect. It was, however, found that median sky re- 
moval was equally well applicable to these observations. See 
Jenness, Lightfoot & Holland (1998) and |Archibald et al. 



(2002|) |br examples of the sky-noise removal algorithm 



3.4 Residual sky removal 

The sky-noise removal algorithm currently assumes that the 
sky signal adds a d.c. offset across the whole array. In some 
extreme cases this is not true and a linear gradient is visible 
across the array. To compensate for this a plane is fitted 
to the regridded image and subtracted before the data are 
processed. For the data considered here this procedure has 
only a minor effect because the source is in the centre of the 
image but is performed for completeness (when the sky-noise 
is a simple offset this step will have no effect). 



regions towards the edge of the map and determining the 
standard deviation of the pixels in the aperture. Bolometers 
with uncharacteristically high noise are often apparent when 
maps are made of faint sources. In order to avoid these and 
other effects like spikes, which are not always easy to remove 
from small datasets, a number of apertures are used and the 
one with the smallest standard deviation is selected to de- 
termine the overall error in the integrated sum (the aperture 
measures the noise-per-pixel). We verified that this proce- 
dure is valid by reducing a number of maps by hand. 

Once the data have been reduced automatically there 
are inevitably a small number of points that must be re- 
moved by hand because of problems with the data acquisi- 
tion, or with the instrument or telescope. The main problems 
are occasional pointing observations where the source is not 
centred on the array, where the dish has been heated by the 
Sun and is grossly out of shape, when SCUBA is not op- 
crating correctly, or when the polarimeter is attached. For 
these cases, the observing logs and final data products are 
inspected and the point is removed if warranted by the cir- 
cumstances. 



3.5 Signal determination 

The penultimate step in the reduction process is to deter- 
mine the signal from the source. In the present case there 
is always a single point-like source at the centre of the map 
and this simplifies the analysis. 

Having first made sure that we have accurately located 
the centroid position there are two ways to determine the 
signal from a source: (1) measure the peak flux, and (2) 
measure the flux in an aperture. Both of these approaches 
are valid, we calculated both, and compare them in Section 

yj 

The first step for case (2) is to select an aperture of 
suitable diameter. For the special case of pointing observa- 
tions the aperture diameter can be no larger than 60 arcsec 
otherwise contamination from the negative images of the 
central source becomes a problem. In general, there are two 
competing considerations, the effects of which have to to 
traded-off to obtain the optimum aperture size for the par- 
ticular problem. On the one hand, it is desirable to use as 
large an aperture as possible so that most of the time de- 
pendent error beam pattern is enclosed for both source and 
calibrator. This is especially important at 450- [im where 
the error lobe contribution typically exceeds 100 per cent 
for chop throws larger than about 2 arcmin. On the other 
hand, increasing the aperture also increases the measured 
noise. The final trade-off is thus between accuracy of calibra- 
tion and signal-to-noise. Experimenting with various aper- 
tures we found that 40 arcsec is a good choice for point-like 
sources. Increasing the aperture to 60 arcsec results in an in- 
crease in signal of only 5-6 per cent at 850- u.m under most 
conditions considered, and correspondingly 13 per cent at 
450- urn whereas the signal-to-noise degrades by a factor of 
2/3 (the error on the integrated sum is proportional to the 
radius of the aperture). Apertures of significantly less than 
40 arcsec produced unacceptable scatter in the signals mea- 
sured for Uranus with the standard deviation increasing by 
almost a factor of 2 (see Table 1) . 

The signal-to-noise is determined by looking at small 



4 RESULTS 



4.1 Flux conversion factor 



The final step in the reduction process is the application of 
a suitable 'flux conversion factor' (FCF) to the reduced data 
so that the measured signal in volts is converted into physi- 
cally useful units such as Janskys. In order to automate the 
reduction of arbitrary data we need to know the FCF history 
of the telescope and instrument. We do this by measuring 
the response of our primary calibrators (Mars and Uranus) 
throughout the period of interest and comparing this with 
the expected flux derived from a standard mode l for plan- 
etary fluxes (|Privctt, Jenness fc Matthews 1998 ; Griffin & 



Orton 1993 



as mentioned in Section 3.5 



Wright 1976|). The FCF can take two forms: 
we can either calibrate on the 



'peak' of a point source giving the traditional radio astron- 
omy units of Jy beam -1 {FCF peak = S peak /V peak where 
Speak is the flux in the beam and V paak is the measured 
peak voltage), or we can 'integrate' in an aperture giving 
units of Jy arcsec -2 (FCF int = Stot/[V in tA] where S to t is 
the total flux of the calibrator, Vi„t is the integrated sum 
in volts and A is the pixel area in arcsec 2 ). In this section 
we discuss the merits of the two schemes by applying the 
automated pipeline reduction to Uranus. 

The results of our Uranus reduction for the 850N and 
450N filters are shown in Figs. |l| and ^| respectively. Simi- 
lar trends are seen for the 450W/850W filter pair, and the 
FCF values are summarized in Table |l|. The top left panel 
of Fig. |l| shows the FCF calculated by summing the signal 
in a 40 arcsec aperture plotted against time for a period of 
2 — 3 yrs. In general the FCF is remarkably stable given 
the variable weather/telescope/instrument parameters over 
this relatively long time period. There is evidence for the 
FCF to be slightly higher during 1997 and early 1998 than 
in subsequent times. It is presently unclear as to why this 
should be the case (if real) but we note that SCUBA was 
warmed up in early 1998 because of varying base tempera- 
ture problems and a similar, small, change in FCF was seen 
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FCF peak / Jy beam 1 V 1 Time of day (HST) / hr 

Figure 1. Flux conversion factors derived for the 850N filter for Uranus. 'FCF int' is the integrated flux conversion factor and 'FCF 
peak' is the flux conversion factor derived from the peak of the source. The vertical dashed lines on the bottom right panel delineate the 
observing shifts of the JCMT. 



following a warm up in the middle of 2000. Also apparent 
from this panel is that a small fraction of the FCF values are 
discrepant by a large amount (> 50 per cent). The effect is 
shown clearly in the top right panel of Fig. ^ where we show 
a histogram of the integrated FCF values. The dotted line 
shows a Gaussian fit to the data giving a standard deviation 
of around 10 per cent as used by Robson et al. (2001). The 
solid line shows the result of fitting a double Gaussian to the 
data (Table ji] lists the average for the filter over the entire 
period) which gives an improved fit. Note the tail of high 
FCF values. 

The lower left panel of Fig. [i] shows the 'integrated' FCF 
versus 'peak' FCF. Variations in the sensitivity of SCUBA 
connected with warm-ups must contribute significantly to 
this general trend shown in this plot. However, on any given 
night we believe that the linear trend, and hence the domi- 
nant source of uncertainty in the calibration of SCUBA data 
is produced by the thermal state of the antenna; it is well 
known that the focus of the antenna changes rapidly around 
sunset and sunrise reflecting its changing shape as it ther- 
mally relaxes. During these periods the observed beam pat- 
tern is also variable because power is removed from the main 
beam and spread out into the error lobes. The fact that we 
see a trend demonstrates that the error beam spreads out 
beyond the 40 arcsec aperture, but the key point to note 
is that whereas the 'peak' FCF varies by about a factor of 
two, the 'integrated' FCF varies by < 50 per cent. Further- 



more, the size of the chop throw can vary the 'peak' FCF 
by as much as 5 — 10 per cent (the beam smears out a small 
amount in the direction of the chop throw and the sensitiv- 
ity is reduced) but affects the 'integrated' FCF much less. 
The best reduction method for our automated approach is 
thus to use the 'integrated' FCF for flux determination. A 
corollary is that photometry mode observations, which are 
made with a single pixel and a 2 arcsec, 3x3 square jiggle 
pattern, and hence are calibrated in an analogous manner 
to the FCF 'peak' method described above, can at times 
be susceptible to large calibration uncertainties. A quanti- 
tative assessment would require detailed analysis of calibra- 
tion data from individual nights which is beyond the scope of 
the automated approach described in this paper, but if high 
values of the FCF are measured then it may well be good ob- 
serving practice to use SCUBA in 'jiggle-map' mode rather 
than photometry mode. Although mapping observations are 
intrinsically less efficient than photometry observations, the 
ability to measure the signal in a suitably sized aperture 
will allow observations to be calibrated accurately during 
the period when the antenna is changing shape. 

The extreme values plotted in the lower left panel of 
Fig. 2 probably occur when the antenna is grossly out of 
shape which may reflect the nature of the observing, for ex- 
ample if observing was carried out in the daytime or if the 
dish was pointed towards the Sun in early evening. These 
data points are calculated from maps in which the beam 
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Figure 2. Same as Fig. 
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but for the 450N filter. 



shape is so obviously poor that they are often removed by 
hand in the post-pipeline phase of the data reduction. Evi- 
dence that this is indeed the case is provided in the bottom 
right panel of Fig. |l| where we plot 'integrated' FCF ver- 
sus time of day Here we see that the majority of high FCF 
values (> 1.3) occur during the first half of the night or af- 
ter sunrise. Note also that the FCF can remain high until 
the middle of the night in some cases but is almost always 
lower than about 1.2 during stable night-time conditions af- 
ter midnight. This may indicate that once the antenna has 
been heated it can take several hours to recover its optimum 
shape. Similar trends are seen at 450-u.m (Fig. ^) where a 
single Gaussian provides a good fit the 'integrated' FCF 
values giving a standard deviation of about 20 per cent. Ex- 
perience shows that the uncertainties derived from our au- 
tomated approach are similar to those applied by observers 
when processing their data 'by hand', demonstrating the va- 
lidity of the method for point-like sources. 



The ratio of 'peak' to 'integrated' FCF provides some 
idea of the magnitude of the error beam within the aper- 
ture at these wavelengths. At 850-u.m the ratio for a 40 arc- 
sec aperture is approximately 250 arcsec 2 beam -1 , an area 
equivalent to a 15 arcsec Gaussian beam (the measured value 
is approximately 14.5 arcsec). At 450- u.m the ratio is approx- 
imately 120 arcsec 2 beam -1 , equivalent to a 10 arcsec beam 
(the main beam is approximately 7.5 arcsec). The average 
error lobe contribution within a 40 arcsec aperture is thus 
about 4 per cent at 850-u.m and 40 per cent at 450- (i.m. 



Table 1. Flux conversion factors determined from Uranus for 
three different apertures. There is evidence for a decrease in FCF 
for 850N in February 1998 (a decrease of approximately 8 per 
cent). This table lists the average for the 850N filter. The Peak 
FCF is only listed once for each filter since the value is not de- 
pendent on aperture. All results are from a combination of chop 
throws between 60 and 120 arcsec. 



Filter Aperture 
diameter 
/ arcsec 



Integrated 
Jy arcsec - 2 V~ 



Peak 
Jy beam -1 V~ 



850N 


20 


1.43 


± 


0.19 


253 


± 


31 




40 


1.00 


± 


0.11 










60 


0.94 


± 


0.10 








850W 


20 


1.19 


± 


0.11 


207 


± 


13 




40 


0.84 


± 


0.04 










60 


0.80 


± 


0.06 








450N 


20 


10.0 


± 


2.6 


855 


± 


260 




40 


6.9 


± 


1.5 










60 


6.1 


± 


1.2 








450W 


20 


4.0 


± 


0.8 


340 


± 


81 




10 


3.2 


± 


0.5 










60 


2.8 


± 


0.5 









4.2 Mars 

In order to test the output from the automated reduction, 
including a test of the validity of our calibration, we tested 
the complete system with observations of Mars. An aperture 
of 60 arcsec was used to maximize signal (signal-to-noise is 
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1.5x10 




2001 



Figure 3. Flux density measured for Mars at 850- (lower curve) 
and 450- ixm (upper curve) . The solid curve indicates the expected 
flux using the model of Wright (1976). Errors are derived by com- 
bining observations taken on the same night. 



Table 2. Coordinates of the calibrator sources as defined in the 
JCMT pointing catalogue. All coordinates are J2000. 





R.A. 


Dec 


HL Tau 


04 31 38.4 


+ 18 13 59. 


CRL618 


04 42 53.60 


+ 36 06 53.7 


OH231. 8+4.2 


07 42 16.94 


- 14 42 49.1 


IRC+10216 


09 47 57.38 


+ 13 16 43.7 


CRL2688 


21 02 18.81 


+ 36 41 37.7 



not an issue) and to guarantee the inclusion of the majority 
of the flux (Mars had a maximum diameter of 15 arcsec dur- 



ing thelperioH) The results are shown in Fig F alo ng with 



the expected flux provided by the fluxes software (Privett 
et al. IS 198) whi ch uses the m odel for mean surface tempera- 
ture derived by Wright (1976) and knowledge of the SCUBA 
filters. The excellent agreement between the observed and 
predicted fluxes provides convincing support for our auto- 
mated reduction procedure and the Mars model. 



4.3 Secondary calibrators 

Now that the system has been shown to generate reason- 
able results we can apply these techniques to archive data. 
Here we present the results from processing of the most com- 
mon JCMT secondary calibrators (see Table ^|); the results 
from extragalactic r adio pointing sources are presented in 
Robson et al. (2001). Table H shows the results of the au- 



tomated reduction of the selected calibrator sources using a 
40 arcsec aperture. These data are selected in the same way 
as for Uranus with only pointing and mapping observations 
retrieved from the data archive. This list does not include 
potential candidates, suc h as VY CMa and TW Hya, men- 
tioned in |3andell (1994] ) simply because the archive does 
not currently contain enough observations to justify reliable 
automated reduction]^] Since all these sources are used for 

2 By the end of the year 2000 there were only 12 suitable obser- 
vations of VY CMa and 20 suitable observations of TW Hya in 



Table 3. Secondary calibrato r fluxes for a 40 arcsec aperture. 
The fluxes in italics are from fcandell (200l| ) and are shown for 
comparison even though they are derived from the peak and will 
therefore not include extended flux. The resu lts for IRC+10216 
and OH231.8 are the average values; see § [4.3] for more details on 
the variability of these sources. Fluxes are determined by fitting 
a Gaussian to the histogram of raw data to determine the peak 
and standard deviation. The number of observations is given in 
brackets. 



850 fim 



450 jim 



CRL618 


4.69 


± 


0.37 (1445) 


12.1 


± 


2 


2 (398) 




4.55 


± 


0.20 




11.5 


± 


1 


5 


CRL2688 


6.39 


± 


0.51 


(223) 


30.9 


± 


3 


8 (57) 




5.9 


± 


0.20 




u 


± 


2 




OH231. 8+4.2 


2.75 


± 


0.44 


(272) 


12.7 


± 


2 


2 (59) 




246 


± 


0.20 




10.5 


± 


1 


2 


HLTau 


2.36 


± 


0.24 


(287) 


9.9 


± 


2 


(83) 




2.3 


± 


0.15 




10.0 


± 


1 


2 


IRC+10216 


8.8 


± 


1.1 


(736) 


17.5 


± 


4 


5 (169) 




6.1 




0.2 




13.0 




1 


8 



pointing there are significantly more observations available 
at 850-um than at 450-um. IRAS 16293-2422 has not been 
included in this analysis beca use of the extended nature of 



the source 
Table I 



Sandell 2001) 



indicates that the automated reduction shows 

especially at 850- 



excellent agreement with Sandell (2001) 
um. The fundamental difference is that the numbers pre- 
sented here are integrated over an aperture and will t herefore 
also include flux from extended emission whereas Sandell 



(2001) calibrated his data in Jy/beam. 



4.3.1 HL Tau 



The f lux density shows excellent agreement with Sandell 
(2001) at both wavelengths indicating that there is no evi- 
dence for extended structure out to 20 arcsec radius when 
performing chopped observations. Fig. |^ shows the light- 
curve. Whilst the plot raises the possibility of source vari- 
ation (possibly with a maximum at 1998.5, amplitude of 
100 mjy and period of 600 d) a chi-square fit indicates that 
the data are consistent with zero variability. 



4.3.2 CRL 618 

CRL 618 is the most popular secondary calibrator in use 
at th e JCMT. The flux density shows excellent agreeement 
with fcandell (200l| ) and Sandell (1991) at both wave lengths 
(the latter taken using UKT14 (Duncan et al. 1990)) indi- 
cating that there is no evidence for extended structure out 
to 20 arcsec radius when perform ing chopped observations. 
Knapp, Sandell & Robson (199J ) indicate a possibility of 



the archive. Of those only 2 of the VY Cma and 6 of the TW Hya 
observations had 450- \aa data. At 850- \aa the results from auto- 
mated reduction were 2.1+0.1 Jy for VY Cma and 1.46+0.2 Jy 
for TW Hya 
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Figure 4. Light curve for HL Tau at 850- um. Errors are derived 
by averaging observations over 5 day periods. The mean flux den- 
sity is indicated by dashed lines. 
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Figure 5. Same as Fig. | but for CRL 618 

source variability in the millimetre region but Fig. [s] shows 
that no obvious variability can be seen at 850-u.m between 
mid-1997 and mid-2001. Apparent short term trends (of the 
order of 0.2 yr) in the data stream are most likely con- 
nected with the observing/reduction methodology or with 
FCF variations due to dish heating rather than to source 
variability. 

4.3.3 OH 231.8+4.2 



We g et excellent agreement with the results of Sandell 



(2001). OH231.8 is a long period variable star and the time 



series show clear evidence for variability at 850- urn. A least- 
squares fit of a sine- wave to the time series indicates a period 
of 630 ±20 d. This is somewhat lower, but not necessar ily in- 
consistent with, that se en in the near-infrared : 650 d (Feast 
et al. 1983;), 684 ± 4 ( |Bowers fc Morris 1984 and 708 ± 6 
flKastner et al. 1992|), but our data are relatively noisy. In- 
deed, [Kastner ct al. (1992 ) suggest that the period may have 



changed slightly over the years and it is therefore possible 
that the period has shortened recently. 

Fig. |^ includes a sine- wave with a period of 630 d of the 
form: 



S{Jy) = 2.75 + 0.4 sin ( 2tt 



Y - 1999.08 
1.75 



(1) 



Figure 6. Light curve for OH231.8 at 850- yxm. Errors are derived 
by averaging observations over 5 day periods. The mean flux den- 
sity is indicated by dashed lines. A sine wave of period 630 days 
is shown. 



S(Jy) = 2.75 + 0.4 sin (2tt 



JD 



2 207 j 



630 



(2) 



where JD is the Julian day. If the sine wave is subtracted 
from the time series data the error in the residual is reduced 
to 0.2 Jy. The uncertainty in the period makes it difficult 
to extrapolate the light curve from earlier studies with any 
accuracy. Therefore it is not possible to reliably determine 
the lag time between the light maximum in the near-infrared 
and that of the submillimetre. The 450-u.m data are too 
noisy to show source variability. 



4.3.4 IRC+ 10216 (CW Leo) 

This is a variable sta r with a period o f 635 d as see n in the 
optical ( lAlksnis 198sh and a t 1.1- mm flSandell 1994| ). When 
compared to Sandell (2001) and Sandell (1994) this source 
is clearly extended in our aperture. Using a source size o f 
8.9 arcsec at 850 -u.m and 4 arc sec at 450-u.m (Sandell 200C) 
the results from Sandell (2001) would be equivalent to 8.4- 
and 16.4-Jy in excellent agreement with our results. 

Fig. [j] shows the measured light-curve at 850-u.m along 
with a sine-wave of period 635 d of the following functional 
form: 



S( Jy) =8.8 + 0.95 sin ( 2tt 



Y - 1999.125 
T74 



where Y is the year, or alternatively: 



S(Jy) =8.8 + 0.95 sin (2tt 



JD 



2451224^ 



635 



(3) 



(4) 



where JD is the Julian day. If the sine-wave is subtracted 
from the time-series data the error in the residual is reduced 
to 0.56 Jy. These data indicate that the source is at a max- 
imum approxima tely 120 d afte r the maximum expected in 
the optical from |Alksnis (1989| ) and 72 d (0.2 yr) a fter the 
maximum expected from the 1.1-mm observations of Sandell 



(1994). It is likely that the discrepancy in phase between the 



where Y is the year, or alternatively: 



850-u.m and 1.1-mm data is not significant given the diffi- 
culty in obtaining the 1.1-mm data with a single pixel in- 
strument. The 450-u.m data are too noisy and do not show 
reliable evidence of variability. 
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UT Date 

Figure 7. Light curve for IRC+10216 at 850-um. Errors are de- 
rived by averaging observations over 5 day periods. The mean 
flux density is indicated by a dashed line. The sine wave has the 
same periodicity as the variation observed in the optical. 



agreement between observed and predicted submillimetre 
fluxes for Mars. It provides a useful tool for investigating 
trends in instrument and telescope parameters - these stud- 
ies are essentially made practical for the first time with our 
automated procedure, and will ultimately lead to a better 
understanding of the uncertainties involved in the calibra- 
tion of data from the JCMT and other submillimetre obser- 
vatories. 

We use our automated procedure to calculate 
light-curves and integrated fluxes for selected, compact 
JCMT secondary calibrators. The periodic variability of 
IRC+10216 is confirmed and there is strong evidence to sup- 
port similar behaviour for OH231.8 at 850- [im; the first time 
this variability has been confirmed in the submillimetre. The 
remaining calibrators have been constant during the period 
1997.5 to 2001.5. 
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U'98 
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Figure 8. Same as Fig. but for CRL 2688 



2001 



4.3.5 CRL 

CRL 2688 h as a 25 per ce nt higher 450- urn flux here than 
presented in [Bandcll (2001 ) suggesting that the source is ex- 
tended at this wavelength. A map was generated combining 
the data from a number of stable nights from the archive to 
test this, taking particular care that any pointing shifts are 
accounted for during the reduction. The measured 'source 
size' of 9.5 arcsec is larger than that measured from CRL 618 
(known to be unres olved) , using the same techniques, of 
8.0 arcsec (see also, Coulson 2001). This corresponds to a 
deconvolved source size of approximately 5 arcsec an d would 
accou nt completely for the lower flux measured by Sandell 
(2001). With a source size of 5 arcsec we would expect a 



much larger discrepancy at 450- \xm than at 850- urn which 
is indeed the case. The time-series in Fig. ^| indicates that 
the flux density has been constant since the start of 1998. 



5 CONCLUSIONS 

This paper presents a means of automating the reduction 
and calibration of 850- and 450- um SCUBA data. The 
limitations of the technique are investigated by analysing 
archival data for the standard submillimetre calibration 
source, Uranus. We show that our method gives excellent 
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